Entanglement generation in persistent current qubits 
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In this paper we investigate the generation of entanglement between two persistent current qubits. 
The qubits are coupled inductively to each other and to a common bias field, which is used to control 
the qubit behaviour and is represented schematically by a linear oscillator mode. We consider the 
use of classical and quantum representations for the qubit control fields and how fiuctuations in the 
control fields tend to suppress entanglement. In particular, we demonstrate how fluctuations in the 
bias fields affect the entanglement generated between persistent current qubits and may limit the 
ability to design practical systems. 

PACS numbers: 03.65.-w, 74.50.+r, 85.25.Dq 



INTRODUCTION 

In this paper, we investigate the use of inductive cou- 
pling to generate entanglement between two persistent 
current qubits. We are particularly interested in the rep- 
resentation of the magnetic bias fields that are used to 
control the behaviour of the qubits, and the requirements 
placed on these fields by the need to generate significant 
levels of entanglement. Given the recent experimental 
results indicating coherent quantum behaviour in super- 
conducting persistent current and other Josephson de- 
vices 0, U 13 1 extension of these systems to arrays 
of coupled qubits for quantum information processing is 
important and timely. Indeed, experiments have already 
been reported using coupled persistent current qubits [J| . 
We demonstrate that the requirements placed on the bias 
fields could present significant obstacles to the use of per- 
sistant current qubits in quantum information process- 
ing. 

We consider two models for the bias field: one clas- 
sical and one quantum mechanical. In each model, the 
field is represented by a lossy linear oscillator, whose res- 
onant frequency and coupling to the qubits can be var- 
ied. When the natural frequency of the field mode is 
significantly lower than the qubit frequencies, the conven- 
tional approach is to treat the bias as a classical variable 
dllb, 4] and to use the expectation value of the screen- 
ing current in the classical equations of motion when the 
bias dynamics are important |^|^. However, where the 
bias field has fluctuations at frequencies that are com- 
parable with the qubit fluctuations, the classical model 
is no longer valid and the quantum mechanical model 
predicts some interesting dynamical behaviour. In addi- 
tion, we derive a set of constraints for the accuracy of the 
bias fields which must be obeyed for a significant amount 
of entanglement to be produced. These constraints may 
limit the entanglement that can be produced in a prac- 
tical system. In particular, we derive restrictions on 




FIG. 1: Schematic diagram of coupled qubit system. 



the coupling between the qubits and the bias fields and 
the operating frequencies of multiple persistent current 
qubits. 



COUPLED PERSISTENT CURRENT QUBITS 

The persistent current qubits studied in this paper 
have been proposed by Orlando et al. ul]. A schematic 
circuit diagram is given in Figure 1, with inductive cou- 
pling between the two qubits and the common control 
field. The qubit inductance is negligible when compared 
with the effective inductance generated by the series 
Josephson junctions in the loop. This means that the 
behaviour of an isolated qubit will tend to be dominated 
by the series Josephson junctions rather than the geomet- 
rical inductance of the ring, Lqu — 10 pH. This allows the 
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circuit to be simplified to a two-state model, correspond- 
ing to current states differing by approximately 600 nA 
Pj. The two-state Hamiltonian for a single qubit, with 
two control fields and ^x2, is given by 0, 



$.2) 



(1) 

where the basis states, {|0), |1)}, are the persistent cur- 
rent states with approximately ±300 nA, is the pri- 
mary bias field for the main ring circuit and is a 
secondary bias field that is used to modulate the critical 
current of the effective Josephson junction formed by the 
two parallel junctions in the smaller secondary ring. The 
matrix elements are given by. 
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and $0 = h/2e — 2 x 10^^^ Wb. The circuit constants 

are taken from p]: ri — 2'kEj^I — — 2r2, si ~ 0, 

ti = 0.001£;j, T] = 3.5, P = 0.8, Ej = 200 GHz, Ec = 
Ej /80. 

Although the energy level separation, and hence the 
dynamics of a single persistent current qubit, is domi- 
nated by the Josephson energy of the junctions in the cir- 
cuit, the inductance is important when determining the 
coupling between the qubit and the external fields and 
between the qubits themselves. For the system shown in 
Figure 1, this gives 
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where Mi is the mutual inductance between the qubits 
and the bias coil and M2 is the mutual inductance be- 
tween the two qubits, $osc is the magnetic flux in the 
shared bias field (which is treated as a linear oscillator 
and characterized by a capacitance Cose and an induc- 
tance Lose)- In the absence of dissipation, the effective 
Hamiltonian for the combined system can be written in 
the form lal, 
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where is an external current used to fix the static 
bias point (the oscillator fluctuates about this point), and 



the coupling coefficients are given by = / LquLosc, 
fii = Ml /Lose and K2 = M|/L^„ = /ij. Each qubit has 
two bias/control fields, ^xii and $3:21 for qubit 1 and 
^xi2 and ^x22 for qubit 2. The primary control fields 
for the qubits, ^xu and $2:12, are common so we put 

^xll = ^x21 = Ml^osc- 

To derive the Hamiltonian for the two qubits, we ex- 
amine the energy of the inductive circuit components. 
The Hamiltonian terms corresponding to the inductive 
energies will have the form. 
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Expanding the second of these expressions, the cross- 
coupling terms between the two qubits has the form: 
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The other terms corre- 
sponding to a shift in the effective self inductance of the 
qubits and cross-coupling between the qubits and the os- 
cillator are subsumed into the F and B terms. The two 
qubit Hamiltonian then has the form, 

Hqui,qu2 *i'£c2i; ^xl2, ^x22) = 



F2 + A12 



- F2 ■ 



■ F2 + ^^12 



(7) 

in the current basis {IO1G2), IO1I2), II1O2), II1I2)} (which 
is used as the computational basis for the purposes 
of this paper), and where F\ = F{^x\\t^x2\\ Bi — 

B{^xll,'^x2l), F2 = F($^i2, $^2), B2 = B{^xl2, ^x22)- 

The A term comes from the qubit-qubit coupling term 
given in equation (6), with A12 = K2LquIqu, where 
lqu — 300 nA is the magnitude of the screening current 
in the qubit logic (persistent current) states. 

We assume that both qubits are identical and we con- 
sider the dynamics of a (common) primary control field 
^ose, keeping the secondary fields <i>a;2i and $3,22 fixed. 
Both the primary and secondary fields are nominally set 
to zero so that the energy eigenstates of the individual 
qubits are symmetric/anti-symmetric superpositions of 
the qubit current states. Initialising the qubits in a cur- 
rent state will produce coherent oscillations at frequen- 
cies around 400 MHz. Although the fields are nominally 
zero, they all include a fixed error and the primary field 
includes the dynamics of the bias circuit. (We use a com- 
mon primary bias for computational simplicity and be- 
cause the decoherence rate will be lower where any noise 
due to the flux bias is the same for each qubit). 

For simplicity, we assume that the couplings are weak, 
typically Ki = 0.002 and ii:2 = 0.01. This means that 
first order coupling terms will be sufficient for most pur- 
poses. Allowing stronger couplings between the qubits 
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and the bias fields could introduce a range of problems: 
difficulties initialising the qubits in a given state since 
the flux and current states are no longer identical (due 
cross-couplings in equation (1)), and the quantum fluc- 
tuations in the bias coil can affect the ability to generate 
entanglement (see below). 



BIAS FIELDS AND DYNAMICS 

We consider two models: one where the control field 
is a noisy classical field and one where it is represented 
by a quantum oscillator. The classical oscillator model 
0, is expected to be valid as long as (i) the typical 
frequency of the oscillator is significantly lower than that 
of qubits, (ii) the (possibly entangled) two-qubit state 
and the oscillator state are separable, and (iii) as long as 
the quantum fluctuations of an equivalent quantum os- 
cillator (approximately given by the width of the energy 
eigenstates in a magnetic flux basis) are small compared 
to the other fluctuations that couple to the qubits. The 
quantum model uses a standard harmonic oscillator basis 
for the control field, and couples via an oscillator flux op- 
erator (formed from the raising and lowering operators, 

and a) in the F and B functions. 

The classical approximation is based on the Born- 
Oppenheimer approximation that is often used in nuclear 
and molecular physics. This removes the dynamics of a 
'fast' degree of freedom by replacing the quantum me- 
chanical operators with their expectation values; thereby 
averaging or integrating out the effect of their dynamics 
on the the other 'slow' degrees of freedom. The details of 
the approximation and the restrictions on its use are more 
fully described in reference 6]. In this case, the qubit be- 
haviour is assumed to be fast compared to the evolution 
of the classical oscillator, and the expectation value of 
the energy is included in the (now classical) Hamiltonian 
given in equation (5). The classical equation of motion 
is then derived in the conventional way using the varia- 
tional derivative with respect to the oscillator magnetic 
flux $osc7 and the energy expectation value becomes the 
expectation value of the combined qubit screening cur- 
rent. Using this approximation, and adding a parallel 
resistance Rose, the equation of motion is given by |^|^, 



+ 



^osc I fl^^^ dt 

= lin +Ml (^ni(/Ul*osc,$2;2l) + (A^I ^osc , $2;22 ) 



(8) 

where the qubit screening currents are calculated from 
the expectation value of the qubit screening current op- 
erators Iqui and /q„2 over the instantaneous wavefunc- 
tion (i.e. a pure state) of the two-qubit state (calcu- 
lated using the time-dependent Schrodinger equation). 
The time-dependent Schrodinger equation is used for the 
qubit evolution in this case because, for simplicity, we 



assume that the dominant source of decohence is the os- 
cillator and any intrinsic dissipation due to emission from 
the qubits in the cavity is comparatively small. However, 
the effect of this emssion process on the behaviour of a 
classical oscillator has been examined elsewhere 0- The 
dissipative term acts as a source of classical fluctuations 
due to Johnson noise in the resistor at finite temperature, 
taken to be T = lOmK which is in line with experimen- 
tal systems. (The noise need not be thermal, but it is a 
useful generic model for experimental noise because elec- 
tronic noise is often characterised in terms of an effective 
noise 'temperature'). 



QUANTUM EVOLUTION AND QUANTUM 
JUMPS 



In the quantum model, the reduced density operator 
for the qubits is estimated using a quantum trajectory 
model: an unravelling of the Markovian Master equation 
that produces individual 'trajectories', which can then 
be averaged over an ensemble to produce an estimate 
of the density operator H i, [13. (A recent compre- 
hensive review of this subject is given in reference 
Each unravelling is equivalent to the Master equation 
when averaged over an ensemble, but corresponds to a 
different measurement interaction at the individual sys- 
tem level |lO| |. For simplicity, we choose the 'quantum 
jumps' model 01 and thermal environment (Lindblad) 
operators for the oscillator described in [l^ . 



Li = [{n + l)uJoscQoscV- 



cQos 



(9) 



where h = [ex.p{hLL>osc/ kT) — is the thermal oscillator 
occupancy, uJosc = 1/ VGoIcLoIc is the resonant frequency 
of the bias feld and the quality factor is given by Qosc — 
uJoscRoscCosc (in this paper we use Qosc = 200 for both 
the classical and quantum models). These environmental 
operators represent emission and absorption of photons 
from the environment by the oscillator mode, which is 
assumed to be the dominant source of dissipation in this 
paper. 

The quantum jump evolution is calculated by numeri- 
cally integrating the full state (describing the qubits and 
the oscillator) over discrete time intervals (of size dt) and 
applying three different evolution operators. In each time 
interval, there is a small (but finite) probability that the 
bias oscillator will emit or absorb a quantum of energy 
from the environment. The probabilities for emission and 
absorption during the time step are found from, 

Pi{dt) = (^L\Li^ dt = [{n+l)uJoscQosc] {a^d) dt (10) 



P2{dt) = {LlL2)dt^[ 



Qo 



osc^osc 
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t)di (11) 
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The jumps are generated stochastically and when a jump 
occurs, a projection operator is applied to the instan- 
taneous state of the system. If an emission occurs, an 
operator Cli{dt) — y/dtLi is applied, lowering the state 
of the oscillator, and if absorption occurs an operator 
Cl2{dt) — \/dtL2 is applied, raising the oscillator state. 
The state is then renormalised. In the absence of a quan- 
tum jump, the evolution of the system is found from the 
non-unitary evolution operator, 

Clo{dt) = l-'-^H-^{L\L,+LlL2) (12) 

The non-unitary term is added to ensure that the evolu- 
tion of the density operator for the coupled system agrees 
with that predicted by the Markovian Master equation, 
when the (pure state) density operators generated from 
an ensemble of individual 'trajectories' are averaged to 
produce an estimate for the mixed state density matrix 
for the whole system ptotai- two qubits and bias oscillator. 
The reduced density operator for the two qubits pi+2 is 
then found by performing a partial trace over the oscilla- 
tor states. (The validity of the quantum jumps approach 
can be checked by selecting a different unravelling for 
the Master equation, and in this paper the results have 
been verified by comparing the behaviour obtained from 
the quantum jumps unravelling to equivalent results ob- 
tained from the quantum state diffusion unravelling 0. 
Quantum state diffusion produces a continuous stochas- 
tic evolution of the quantum state, and can be shown 
to correspond to a unit-efRciency heterodyne detection 
measurement process [l3)- 

RESULTS 

In each model, quantum and classical, the initial con- 
ditions set for the oscillator are a thermalised state: a 
thermal quantum state for the quantum oscillator and 
an initial condition generated from the classical equation 
of motion, which has been allowed to come into equi- 
librium with the thermal noise by numerically integrat- 
ing its behaviour prior to initialisation using a different 
realisation of the noise for each qubit trajectory calcu- 
lation. The qubits are initialised in a product of pure 
current states. The initial states of the qubits are cho- 
sen to be either in-phase (IO1O2)) or anti-phase (|0il2)). 
That is, the coherent oscillations induced by initialising 
each qubit in a current state are initially in phase with 
each other or in anti-phase. Each individual quantum 
trajectory is calculated using the same initial conditions 
for the qubit states with different static bias errors for the 
control fields. The static bias errors are fixed for each tra- 
jectory and represent the accuracy with which the fields 
might be set in an experiment. The size of these static 
bias errors is found to be crucial to the generation of 
useable entangement between the qubits. The entangle- 



ment is characterised in terms of the concurrence, which 
is widely used in quantum information processing for bi- 
partite systems !l3|. The concurrence for the two-qubit 
mixed state density matrix pi+2 is defined as 

C{pi+2) = max{0, v/a7 - - ^/X^ - \/%:} 

where \/A]", . . . , ^/Xa are the eigenvalues of the matrix 
Pi+2{o'y (8) o'y) Pi+2i^y ® ^v) nondecreasing order and 
CTy is a Pauli spin matrix il4j . Although the concurrence 
is used here, other measures of entanglement can be cal- 
culated from this: e.g. the entanglement of formation 
Ep can be calculated from 



Ef{p,+2) - h (^-[1 + - C(pi+2)2] j 

where 

h{x) — ~x log2 X — (1 — x) log2(l — x) 

at least for this two-qubit system. The mixedness of the 
resultant two-qubit state is given by the von Neumann 
entropy of the reduced density operator Q: S{pi+2) = 
Tr{pi+2 log4 Pi+2)- The logarithm is taken to base 4 be- 
cause the qubit states exist in a four-dimensional Hilbert 
space, giving a mixedness parameter that varies between 
zero and one. 

Figure 2 shows the concurrence of the two-qubit mixed 
states for bias errors of A^^: = 10~^$o (1 cr, Gaussian) 
for both models for an oscillator frequency of 100 MHz. 
We see that the concurrence oscillates and decays grad- 
ually in time, mainly due to dephasing between the two 
qubits originating from the static error in the bias fields, 
and approaches a maximially entangled mixed state in 
each oscillation. By varying the size of the errors, we 
find that the concurrence is quite sensitive to errors in 
the static fields, bias errors larger than about 5 x 10~^$o 
do not lead to any useful entanglement in the qubits - the 
concurrence is less than 0.3 at all times and decays very 
rapidly. (Other calculations away from the minimum 
splitting point of the qubits, where the frequency dif- 
ferences between energy eigenstates are higher, indicate 
that the requirements on the bias fields are even more 
demanding. Where the qubit frequencies are around 1-2 
GHz, useful entanglement is only generated if the bias 
errors are less than about ^ 10~^3>o)- 

By increasing the couplings between the two qubits, 
the rate at which the two qubits become entangled can be 
increased. However, this could lead to problems when ini- 
tialising the qubit states since the qubit current and flux 
are mixed by the inductive coupling. Care is required to 
ensure that the initialisation process projects the states 
onto the correct basis. (Slight differences in the initial 
states will affect the entanglement in the mixed state, 
but this is not explicitly considered here) . Increasing the 
coupling between the oscillator and the qubits is likely to 
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(a) 1 




FIG. 2: (a) Concurrence versus number of oscillator cycles 
(lOOMHz) for both the classical (blue) and quantum (green) 
field models, (b) Concurrence versus von Neumann entropy 
for mixed state (blue) corresponding to the quantum model 
shown in (a), with two types of maximally entangled mixed 
states (Rank 3 (green) and Werner states (red) (iJI)- 



lead to additional problems. Although the classical os- 
cillator model contains noise due to thermal fluctuations 
and dissipation from the finite quality factor of the oscil- 
lator, a quantum oscillator also includes quantum fluctu- 
ations. The differences in Figure 2(a) for the quantum 
and classical oscillator models are due to the compara- 
tively low Q value used and the quantum fluctuations 
coupling across to the qubits. Increasing the oscillator 
Q and/or reducing the oscillator frequency improves the 
agreement between classical and quantum models. The 
coupling between the bias and the qubits is sufficiently 
small for the entanglement between the oscillator and 
qubits to be negligible. 

This raises an interesting point: what happens when 
the quantum fluctuations in the oscillator coupling across 
to the qubits are comparable with the static bias errors? 
The size of the quantum fluctuations in the oscillator 
can be estimated from the flux width of the harmonic 




Number of Oscillator Cycles 



FIG. 3: Concurrence versus number of oscillator cycles 
(lOOMHz) for quantum field model with different coupling 
strengths in the absence of static bias errors: A'l = 0.002 
(blue), Ki_ = 0.005 (green), Kx = 0.01 (red). 

oscillator states. Using the width of the oscillator states 
and the coupling coefficient, the approximate size of the 
fluctuations that couple to the qubits will be, 

//iA$o,, ~ Kiy/2huJoscLqu 1.2 X 10-6$o (13) 

for Ki = 0.002 and an oscillator frequency of 100 MHz. 
The size of the quantum fluctuations that couple across 
varies linearly with Ki and as the square root of the 
frequency. This means that the field frequencies must 
be very low if strong field-qubit couplings are to be 
used. Keeping the frequency constant and increasing 
the oscillator-qubit coupling we find that as soon as 
the quantum fluctuations become comparable with the 
constraints on the bias errors, entanglement is effec- 
tively lost (~ 10~^$o for the cases considered here and 
^ 10~^$o for qubits biased away from the minimum 
splitting point). Even in the absence of the static bias 
errors, the quantum noise will affect the generation of en- 
tangement between the two qubits. Figure 3 shows the 
effect of increasing the size of the quantum noise by in- 
creasing the coupling to the bias field. For a 100 MHz 
oscillator and a coupling of Ki = 0.01, giving fiuctua- 
tions /iiA$osc — 6 X 10~^$o, the entanglement between 
the two qubits is lost very quickly. That this is due to 
quantum fluctuations, rather than the change in coupling 
strength alone, can be verified by simultaneously chang- 
ing the oscillator frequency and the coupling strength 
keeping the size of the quantum noise given by equation 
(13) fixed. 

The quantum fiuctuations effectively limit the oper- 
ating frequency of persistent current qubits as quantum 
processing devices, because the operating frequency must 
be lower than the frequency at which the bias fields may 
be manipulated, which is determined by the frequency 
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Number of Oscillator Cycles 60 




FIG. 4: (a) Concurrence versus number of oscillator cycles 
(380 MHz) using the quantum model for qubits initialised in 
the IO1I2) state (blue) and the IO1O2) state (green); (b) Con- 
currence versus von Neumann entropy for qubits initialised in 
the IO1O2) state (blue), with two types of maximally entangled 
mixed states (Rank 3 (green) and Werner states (red) [ij); 
(c) as (b) for the ]0il2) state. 

and the quality factor of the bias circuit. Increasing the 
operating frequency of the device, and keeping the fluctu- 
ations below the required level, would mean reducing the 
coupling between the qubits and the applied field, which 
might make it difficult to address individual elements of 
an array of qubits. 

In spite of the possible difficulties in biasing and ad- 
dressing individual qubits within an array, there are some 
aspects of the behaviour of this tripartite system that are 
worth investigating further. In particular, in situations 
where the frequencies of the oscillator and the qubits 
are not widely separated, it is possible to generate in- 
teresting evolution whereby the currents flowing in the 
qubits excite oscillations in the bias field, which modifies 
the behaviour of the qubits. Figure 4 shows examples 
for a quantum oscillator with a natural frequency of 380 
MHz, for both the in-phase and anti-phase initial states 
(all other parameters are identical to Figure 2). For the 
IO1O2) initial state, the screening currents flowing in the 
qubits add in-phase. The net current coupled to the oscil- 
lator acts as a sinusoidal drive which excites oscillations 
in the bias field. In the IO1I2) initial state the net cur- 
rent coupling to the oscillator is close to zero initially and 
the concurrence oscillations are far more regular. There 
are significant differences in the concurrence and the von 



Neumann entanglement between the two cases. In par- 
ticular, the entanglement persists for longer for the IO1I2) 
initial state and the von Neumann entropy exhibits some 
large scale oscillations, as shown by the 'loops' in Figures 
4(b) and 4(c). These oscillations in entropy correspond 
to points where the oscillations in the bias field are at 
their largest. The oscillations in the field shift the bias 
point of both of the qubits, which accentuates the natural 
dephasing between the qubits. As the relative phase of 
the qubit oscillations changes, the net current coupling to 
the oscillator changes, and - as they approach the anti- 
phase state - the net current coupling to the oscillator 
falls and the oscillations in the field reduce, thereby sta- 
bilising the relative phase of the coherent oscillations. Al- 
though these phase slips occur for the anti-phase initial 
state, they are more evident for the in-phase initial state 
and are responsible for the rapid decay of entanglement 
in this case and are the dominant source of decoherence. 



CONCLUSIONS 

In this paper, we have discussed a coupled system con- 
sisting of two persistent current qubits and a linear os- 
cillator, representing one of the qubit control fields. We 
have examined the generation of entanglement between 
the two qubits in the presence of a dynamical bias field, 
and have shown that classical models are approximately 
valid for low frequency fields with high quality factors. 
We have used the models to set constraints on the accu- 
racy of the applied control fields. The static bias errors 
must be less than about 5 x 10~^$o for the parameters 
used in this paper. We have also considered cases where 
the underlying quantum fluctuations in the applied field 
are significant, and have used this to derive a constraint 
that relates the coupling between the bias field and the 
qubits and the frequencies present in the bias fields. If 
these constraints are not met, the useful entanglement 
between the two qubits is effectively lost. This could af- 
fect the use of these devices in a practical quantum pro- 
cessing system, placing severe demands on the accuracy 
of the static control fields and limiting the operating fre- 
quencies of these devices. However, we have found that, 
where the frequencies of the applied fields and the qubits 
are comparable, some interesting dynamical behaviour 
can be produced by the back reaction of the qubits on 
the applied field. 
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